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We make an extensive use of BLM approach to study details of predicting higher order corrections 
based on the approach. This way we are able to test two procedures to improve the prediction 
process. Beside the main line of the two procedures it is found out that overall normalization could 
change BLM patterns effectively. Finally we try to find out whether a BLM pattern is sufBcient for 
a prediction or not, and how one should use such a pattern. 
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I. INTRODUCTION 

Optimization procedures provide a good framework in 
investigating the behavior of perturbative descriptions of 
gauge field theories' observables. Perturbative analysis 
of these theories encounters the problem of scheme-scale 
ambiguity. An appropriate approach to this problem 
should give warnings in cases where it should not be ap- 
plied, or we should have a mechanism to check the valid- 
ity. BLM approach [l[ contains such a mechanism, it is 
not applicable in analyzing processes containing gluon- 
gluon interactions in leading order since we should not 
absorb all the flavor dependent part into the running cou- 
pling constant in these cases. If we use the approach to 
resolve scheme ambiguity and find physical schemes, such 
processes will automatically be eliminated. Strange per- 
turbative series would be the outcome of using them as 
the physical scheme. The mathematical form of pertur- 
bative series suggests a kind of symmetry and relation be- 
tween observables but without considering optimization 
processes the real hidden relation can not be revealed. 
BLM approach is particularly able to reveal commensu- 
rate scale relations between observables. This flexi- 
bility of BLM approach can provide us with very simple 
relations such as the generalized Crewther relation ^] 
which can be used to test some serious aspects of QFT 
such as the violation of conformal symmetry due to the 
renormalization procedure. 

One can use optimization procedures to predict higher- 
order coefficients. Clearly any kind of prediction strongly 
assumes existence of patterns. A pattern at least must 
not have mathematical defects. However the quality of 
a proposed pattern is usually determined through the 
comparison of the predicted and exact coefficients. For 
a flexible approach this comparison could be a guideline 
in improving the approach. As a matter of fact what 
should be followed here rather than improving an ap- 
proach is finding patterns. Fortunately we do not have 
many choices for the first step. The first step towards 
these aims is to make the prediction rebuild exact values. 
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This way we encounter a new face of the observable. This 
mathematical adaptation can be interpreted as isolating 
parts of the observable. These isolated parts carry the 
physical considerations which supports the optimization 
procedure as parts of them. The set of isolated parts 
can be called a pattern. However this is not the only 
way to find patterns. Regarding the approach extension 
and improvement, we introduce biased BLM approaches. 
Apparently a normal BLM approach cannot correctly an- 
alyze gluon-gluon interactions in leading order. Here the 
suggestion of finding BLM generated parts (the isolated 
parts generated through BLM prediction comparison) in- 
side the observable and using the new obtained parame- 
ters instead of the original ones constitutes an immediate 
biasing. Also we introduce a second-class biasing where 
the leading-order is totally dismissed from prediction. 

From a more illuminative point of view, we check the 
potential of the approach in revealing a pattern for the 
observable. In the case of BLM this has been preferred to 
finding a way to separate gluon-gluon interactions from 
the vacuum polarizations. If the approach can be ad- 
justed to find a satisfying pattern, it has jumped over 
the obstacle. 

Following the route of finding patterns will lead us 
through a pile of equations to a simple sufficient con- 
dition for the relations that connects observables. The 
condition is a sufficient one if someone wants to get new 
results from prediction through a second observable. 

For the convenience of the reader: firstly in this way 
of finding patterns and biasing we deal with systems of 
equations that are built upon comparisons. These com- 
parisons may aim different types of parameters in the pat- 
tern such as the perturbative series coefficients or solely 
the local symmetry parameters. Secondly we make no- 
tices anywhere one may find himself in a state to perform 
a selection of equations or parameters such as when deal- 
ing with unequal number of equations and parameters. 

In the beginning of Sec|lT]we will discuss primary as- 
pects of prediction. Deviation pattern approach is intro- 
duced in Sec lII Al accompanied by its details and uses for 
Higgs boson decay widths, Bjorken sum rule and Adler 
function and also some complementary discussions on 
prediction. In this section we also consider patterns de- 
veloped from the combination of different patterns. As 
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we go further we will try take advantage of a second ob- 
servable for finding patterns in Sec lIIBI and introduce 
indirect relations. 



II. THE PROCEDURES 

We consider a RG-invariant perturbative expansion 
i?(g2) = ro + ria,(g2) + ^^or,a/(Q2) where a, = ^ 

(cks is the RG-invariant effective coupling constant) and 

= Sfc=o ^ik'iT'f for any effective charge or MS- like 
scheme. Predicting any coefficient using BLM ap- 
proach results in a series of number of flavors n f for that 
order that depends on all the lower order coefficients. For 
instance the predicted third order coefficient is 

^_121 C|r|i _ 17C|^ _ 11 Ca ^20^21 
16 r|, n 8 Tf^^^ 2 Tf ri 

+ (2^^ + -r,,CA + -r2iCF)nf + '^n). (1) 
ri 4 4 ri ■' 

Here we have used the 4-dimensional /3-function of the 
massless MS-like schemes. The nth order predicted co- 
efficient r^*^ — X]fc=o '''^k^^f appears when re-expanding 
the BLM scale-fixed [n — l)th order series 

i? = ro + ria(g*') -|- f2a'(Q*') + . . . + 7v,_ia"-i(g*'). 

g* is to absorb all vacuum polarization. Re-expanding 
the above series happens through the nth order running 
of the couplant: 

a(g*') = a(g2) - a2(g2)/3o£ + a3(g2)(/32^2 _ ^^^^ 

+a4(g2)(-/33^' + ^PoPif - M + ... + 0(a"+i), 

where £ = ln{Q*^/Q^) and it may occur in multiple steps 
with appropriate is for a multi-scale version of the BLM 
scale-fixed series. We will discuss more details of BLM 
extensions later. 

One more preliminary review on overall normaliza- 
tion and Casimir operators is needed before going on to 
the procedures. Local gauge invariance constructs Yang- 
Mills theories i.e., non-Abelian gauge theories. Such lo- 
cal symmetries suites well to compact semi-simple Lie 
algebras. For such an algebra one can find the overall 
normalization that can determine the structure of 
the algebra and its representations. The tensor prod- 
uct of generators of the color SU{N)s of two particles, 
T^T^ , gives the taste of color interaction between them. 
This operator is surprisingly an invariant of the tensor 
product space. The Casimir operator T^, where T^s are 
generators of the product space, related to this inter- 
action is a number depending on the dimension of the 
representation. Ca and Cf are quadratic Casimir op- 
erators of the adjoint and fundamental representations 
transforming gauge and fermion fields. The representa- 
tion's generators scale factor i.e., the trace normalization 



of the representation is Tf for the fundamental repre- 
sentation. Obviously setting Na, the number of gener- 
ators of the group, Cf, Tf, and d'^'^'^, the higher order 
group invariants, to 1 and the parameters of the adjoint 
representation to will separate the U{1) factor of the 
group. The dynamics of the quenched approximation of 
the non-Abelian part is revealed in the limit nf = for 
all physical energies while the conformal invariant part- 
ner of the symmetry group sits where the /3-function is 
0. 



A. Deviation 

If BLM prediction pattern is some good pattern for 
the perturbative series, one can rewrite the whole se- 
ries in terms of just ri and r2 coefficients. The quality 
of the pattern will depend on these parameters beside 
the BLM approach itself or truthfully the mathematical 
mechanism employed for it. Anyway if BLM approach is 
able to discover the pattern, it will not be so surprising 
if the original ri and r2 do not be the ones that develop 
the pattern. A straightforward step would be to bor- 
row the ri and r2 that form the real since it is the 
most important part of the pattern. Yet the procedure 
would be flexible enough to be able to handle exceptional 
but important cases such as the ones with large higher- 
order corrections. We could discuss the situation also as 
if the weights of ri and r2 parts should be other values 
in Eq. ([T]) or other combinations of weights of lower-order 
coefficients are appropriate. 

In general a BLM pattern contains several coupled 
equations corresponding to all we know about the per- 
turbative series. These equations are acquired by com- 
paring the predicted coefficient and the exact ones. For 
instance when one knows in addition to and wants 
to adapt the prediction pattern so that all the details are 
involved in it, he will notice that r^^"^ depends on all the 6 
parts of lower order coefficients {ri, r2o, r2i,r3o, ^31, r32}. 
But the 4th order comparison(r4'^'^ = r^) yields 4 equa- 
tions. Considering the 3rd order comparison, we see that 
7 equations are controlled by just 6 variables. This in- 
dicates that one of these 7 parts is generated automat- 
ically through the adaptation of the other parts. The 
more information, the more automatic generation. So 
here one is forced to do a selection of equations consid- 
ering the importance of the corresponding comparison. 
For instance one may find out that r^^'^ plays the least 
important role in the deviation of r'^^'^ from r4, the com- 
parison r^g*^ = will be easily removed from the list 
of equations. The same works for other selections of pa- 
rameters too. Deviation Pattern procedure obtains new 
parameters through a selection of comparisons and then 
replaces the original parameters with the new ones to 
build the higher-order coefficients. 

It is possible to involve Adler transformation into the 
above procedure by performing the prediction or compar- 
ison for the corresponding Adler function. Some notes 
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regrading perturbative aspects of Adler function will fol- 
low. The real function R{s) where s denotes the center- 
of-mass energy squared is related to the Adler function 
D(Q2) through 



ds 



The reverse transformation rebuilds R{s) as follows: 



R{s) 



2^ 



s+ze 



Dpti-z). 



One may use the RG-invariant perturbative approxima- 
tion of the Adler D-function Dpt{Q'^) = X^n dnd^iQ^) in 
the reverse transformation. This will generate the time- 
like RG-invariant effective coupling oi{s) and a set of pip- 
izated functions. Certainly we deal with the analytical 
continuations [t'I of Euclidean perturbative expansions 
when using the reverse transformation. This results in 
pipizated expressions [1, Q . The non-power sets of pip- 
izated functions and analyticized functions replace time- 
like and space-like couplings which regenerate the per- 
turbative series in terms of a self-consistent scheme. 



1. Higgs decay widths 

As an important and interesting illustration of the De- 
viation Pattern procedure (DPA), we first consider Higgs 
boson decays [10|] into bottom quarks [HI, [11] and gluons 
[isj . The hadronic decay width of the Higgs boson seems 
to have large QCD corrections to the Born approxima- 



tion. Considering Mh >> 2mq results in 



T{H ^ qq) 



for the decay width, where the term to the left of R 
is ^Born{H — > qq). The absorptive part of the corre- 
sponding two point correlator R generates the QCD cor- 
rections. Gluonic decay width takes the following form 
based on the fact that top quarks plays the main role in 
coupling of gluons to Higgs boson. 



99) 



^^^-Ci2lmH^^(q2), 



M 



H 



where Ci carries the top mass dependence and H'^'^ is the 
induced vacuum polarization by the renormalized gluon 
operator. The gluonic decay width is factorized by the 
if-factor, T{H gg) = KTBomiH gg). 

Now we are in a position that we should consider gener- 
alizations of the BLM approach. For a historical overview 
of the extensions to the approach see [ll] and references 
therein. Our emphasis here is on the single-scale exten- 
sion [ll] used in and the multi-scale one developed in 
[2|. Scale or scales are in need of n{n — l)/2 parameters 
to be able to absorb all vacuum polarization insertions in 
the nth order truncated series. The single-scale extension 
uses these parameters at once. However the multi-scale 
version takes them step by step with a specific theoreti- 
cal consideration. Categorizing corrections forms the ba- 
sis of some of the other generalizations such as IJ] also. 
Taking advantage of the single-scale I results in 



pre _ 187£m. ^03_C| 51 C\Cf 121C|r2o 3993 C|r2i llCA^rw 
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+ [i^CA + -CF)r32 - {^CaTf + -CFTF)r2i - {-Ca - ^^f)— + 2—— - 3— p + 2^:^]^^ 

+ (-4 + 2^W, (2) 



as the predicted 4th order coefficient. Strange thing is 
the appearance of Tf in the denominator of Eq. ([T]) and 
Eq.(l2]). This is even more unusual for the multi-scale 



based prediction where the common factor rl{4TFr 20 + 
llCAr2i)Tp sits as the denominator of the whole r^'^- 
Apparently in both cases r^^*^ does not posses a nor- 
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mal coefScient structure. A normal structure is con- 
structed through a combination of /3-coefRcients muhi- 
phcations. For re-expansion we are using a version of ^ 
which contains n[n — l)/2 parameters for the nth order 
series. One suggestion might be to add n extra param- 
eters to t or ^s when re-expanding to the {n + l)th or- 
der series and see whether a standard structure of 
could determine these parameters. Perturbative series 
parameters {ri,r2,r3} presented in Eq.((T]) and Eq.® 
intrinsically contains local symmetry group parameters 
{CaiCp,Tp,...}. So before following the above sug- 
gestion, relations must be rewritten in terms of new se- 
ries parameters which are group parameters independent. 
However we will get around this and take advantage of 
the single-scale extension in the following predictions. 

In this way the normal BLM prediction generates the 
following series as the 4th order coefficient of R compar- 
ing to the calculated one: 

ff'== 1103.4 - 268.480n/ -I- IQ.Qln) - 0.046^, 
f4 = 39.354 - 220.943n/ -I- 9.685nJ - 0.020n^. 

If we perform comparison and prediction for the corre- 
sponding Adler function, we will have: 



-pre 



-1627.6 + 191.328?!/ - 11.41?!^ + 0.245n|^, 



which indicates an improvement. 

The normal prediction for the 4th order coefficient of 
the K factor reads: 



It is more evident in the /3-representation of the coeffi- 
cients: 



r^-,= 81.04- 



259.09;3o + 25.09;3^ 
n = 120.0 - 261. 26/3o + 75. 95/3o 



4.03/32 



U.mpl -f 4.03/32. 



Thegenuine /3-representation and seBLM are developed 
in [1^ but here we have used a different combination of 
/3- coefficients for the /3-representation. As we go further 
and employ DPA, normal prediction will be closer to the 
calculated result. So the first term in its /3-representation 
is near r4's. At the same time the biased prediction de- 
creases its slope and reaches —54.5961 for = 6 in com- 
parison to —382.9029 of the last step while the calculated 
value is —37.2079. Taking advantage of the Adler func- 
tion makes the biased prediction get away from but 
the normal prediction continues decreasing its slope and 
this time gets really close to as n/ increases: 

rf" = 86.31 - 112.75;3o + 29.87;3^ - 2.49/3^ -|- 4.61/32, 

which reads —35.2114 for rt/ = 6. 

If we establish the ratio as(Q^)^, normal prediction 
will generate a result diverging from the calculated coeffi- 
cient as n/ increases. As we use DPA it suddenly changes 
its orientation towards the calculated result. Utilizing 
DPA with the Adler transformation will make them be- 
have more similarly. 



^= 2219.878 - 901.698n/ + 53.524n^ - 0.516n), 
33 = 3372.073 - 866.588n/ -I- 48.088nJ - 0.538n^. 

Performing comparison and prediction for the corre- 
sponding Adler function improves the result: 

gl''" = 3783.694 - 792.353n/ -I- 42.653nJ - 0.581^, 

whereas preventing the first term of the perturbative se- 
ries affect the prediction during re-expansion (a biased 
BLM prediction [l6|) alongside the above procedure re- 
sults in a better behavior for higher number of flavors: 



gf" = 533.831 



189.814n/ -f 15.738n^ - 0.318n^. 



Biased prediction refers to this kind of biasing. 

Briefly normal BLM prediction for both the decays re- 
sults in patterns whose deviations from the calculated 
coefficients do not have serious nf dependences. Spe- 
cially for R it is so clear that the major part of deviation 
is due to the first part of 4th order coefficient. In the 
case of K factor we have two large parts. The deviation 
is shared among both of them, however again the first 
part plays the main role. DPA will make r^'^ recover f4 
for higher values of n/ , whilst this happens for g^'^ when 
DPA is utilized by a second-class biasing. 

If we consider the ratio CsfQ^)-^, a normal BLM pre- 
diction will indicate the same behavior as before. This is 
whilst a biased one shows a tendency between the pre- 
dicted and calculated coefficients for higher values of n/ . 



2. Adler function and Bjorken sum rule 

The vacuum polarization induced by vector current 
— tpi"fn'4'i is the time-ordered correlation function 

(q^q. - q^9,.MQ^) = i J d^xe'''--{0 \ T[j^{x)MO)] I 0), 

where — — g^. The Adler function corresponding 
to the current takes form through the scalar correlator 

n(Q2), 

DiQ^) = -127r2Q2_l_n(Q2)_ 

Integration over polarized proton and neutron structure 
functions g^'"' is related to the Bjorken polarized sum 
rule quantity dsj through 

rr"- [' dxislix, Q')~gl^{x,Q')) 
Jo 

= - '^BjiQ )) + 2^ Q2j-^2 ' 

where (1 — rfsj) denotes the coefficient function Cbjp 
and gA is the charge corresponding to the axial vector 
current of the nucleon. Recently the coefiicient func- 
tion C^^^{Q^) and the non-singlet component of the 
Adler function D^^{Q'^) [l3| have been calculated to 
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order a^. They are based on the constraints due to 
the special form of the generalized Crewther relation 
D{as)C''^Pias) = 1 + ^EiKiai [H [H (/3(a,) is 
the gauge group /3-function) and the color structure of 
the coefficients of these functions [l^ . It is worth notify- 
ing that more recently the generalized Crewther relation 
(CR) has been expressed in terms of new /^-independent 
polynomials by factorizing multiple powers of /3-function 



D{as)C''^P{as) ^l + Y,(P^ \ VM. (3) 

n>l ^ / 

Adapting the old and new structures of CR reveals the 
structure of Vn coefficients. 

Prediction results for C^^f and D'^^ are summarized 
in Table m For C^^^ it seems that a major part of devi- 



TABLE I. C'^^P and L>^^ prediction results. 
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7.69 


0.10 













DPA. 

DPA utilized by Adlcr function. 



ation for the normal prediction is due to the first term. 
But unlike the case of Higgs decay into bottom quarks the 
gap between prediction and calculation is not so large. It 
is whilst D^^ is completely different. The prediction and 
calculated results may seem really close to each other but 
their /3-representations are really far. There seems to be 
a similarity between D^^'s and Higgs decay into gluons 
width's deviation patterns, both of them have extremum 
near Uf = 3. DPA makes an improvement in C^^^ pre- 
diction for higher values of and provides a slight im- 
provement in the /3-representation. A good improvement 
in the first term of /3-representation for D'^^ is visible 
in its DPA prediction. At last taking advantage of Adler 
transformation will not improve C^^^ prediction. 

The results of predictions for the ratio fjBjp^2-^ could 
be illustrated in Table |TT1 The deviation pattern of the 



TABLE II. ^BjF prediction results. 
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1925.9 


-463.57 


24.11 


-0.2872 


239.47 




1409.6 


-332.47 


16.05 


-0.1302 


397.27 


dcY^^ 


111.29 


-140.82 


7.128 


-0.0043 


397.27 


del'"'''* 


2444.4 


-572.69 


31.02 


-0.3415 


595.05 



ratio does not contain any extremum and is decreasing for 
higher values of n/. DPA prediction is good for n/ = 6 
and making use of Adler transformation will not help. It 



seems that C^^P in this ratio affects parts of D^^ for 
which the D^^ pattern was weakened while employing 
DPA. It makes the new predictions better for ti/ = 6. 

Table Hn] shows the results for the ratio ttn^t^iy- 



TABLE III. prediction results. 
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cd'r 


-54.66 


33.49 


-2.26 


0.0319 


-50.64 


cdY^^ 


-101.41 


30.14 


-1.67 


0.0153 


-66.90 


cdl"^ 


-37.55 


18.99 


-1.02 


0.0028 


-66.90 


cdT"^''^ 


-172.83 


45.62 


-3.03 


0.0379 


-150.66 



Existence of D ° in this ratio absorbs the converging 
behavior of DPA prediction for higher values of n j . This 
converts it to a general improvement along all values of 
rij. The corresponding results for the biased BLM pre- 
diction are tabulated in Table IIVI This indicates real 



TABLE IV. biased prediction results 
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52.07 


13.64 


-1.40 
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-72.32 


cdl"^ 
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-21.23 


0.41 





-140.29 
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461.62 


-73.41 


3.19 


-0.0409 


-140.29 




-172.83 


45.62 


-3.03 


0.0379 


-150.66 



improvements in the /3-representation. 



3. Quadratic Casimir invariants 

Each part in the first term of a normal BLM pattern 
Eq. ([2]) contains at least one factor Ca ■ So for the case of 
Higgs decay into bottom quarks setting Ca = improves 
the normal prediction significantly along all values of 
specifically for lower ones. 

Working on the ratio as{Q'^)^ and using {?'i,r2} that 
came out of DPA, one will find out the sensitivity of the 
rij part of the 4th order prediction to Ca in such a way 
that setting Ca — 0, Cp — 1, and Tp — 1 results in a 
0.04 relative error in predicting this part. 

For the ratio as(Q^)-^ an interesting improvement oc- 
curs when choosing the canonical i.e., conventional form 
of the quadratic Casimir invariants for SU{2). More 
freedom is obtained by departing from the conventional 
choice to the overall normalization p2| | where parame- 
ters b and N take control of the normalization. This also 
indicates that SU{2) is a better choice for this ratio. 

Adapting the first part of the 4th order coefficient of 
(jBjp using the overall normalization indicates that 
SU{2) improves the part and SU{A) improves the nj: 
part in comparison to SU{3) while SU{3) still provides 
a better combination. 
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A very exciting result is obtained when performing the 
last procedure for D^^ . In this case SU{2) gives us such 
a combination of and nj parts that makes the BLM 
pattern fit to the calculated terms in the whole region of 
nj. It is in spite of the slight separation occurred in all 
the previous cases for higher values of n/. 

One should pay attention to these parameters of gauge 
symmetry because their effects are noticeable on the 
problem and specially on an approach which seems to 
be related to the concept of conformal symmetry. 



B. Second observable 

In the previous section, we tried to find parts of the 
observable that generate a BLM pattern. It is possible 
to write the observable in a complete BLM pattern form 
in terms of an appropriate effective charge. This means 
an absolute freedom in choosing the first two coefficients 
ei, 62 or equivalently bi, 62: 



R = eia2 
02= bia + b2a? + 630"^ 



4 

b^a^ H 



As in Eas. (|ll2l) we will have the following relations for 
coefficients: 



__121Cie|i _ 17 Cl 
16 n ei 8 Tf''^^ 



11 Ca 620621 
2 Tf 61 



(2 



£20^21 



61 4 
. 11 

64 = ^620631 + 

6iJf 2 

Ci 627e|i _ 187 
T| ^ 32 61 

17 



e2iCA + -;e2iCF)nf + 



61 



33 33 
— 6206210^ — —r^30S2l\ 
o 4 



Cl 203 

^632] + ^Yl52 



621 



-631] 



Cl ,1089 



6iT/^ 64 



363 



621631] + 



(4) 



so there is an infinite number of observables in terms of 
which a BLM pattern is obtained for R. 

A very first choice for ei, 62 would be ri, r2. For R the 
predicted result suffers the same problem as the normal 
BLM prediction i.e., a large deviation due to the first 
term. An impossible state would be the generation of 
both 6 and b through BLM pattern. 

Obviously we cannot use strict constraints to obtain 
some 61, &2 because these cannot be partners to the first 
constraint due to the BLM pattern. For instance if 02 
is determined as the Adler function corresponding to i?, 
61 will not be able to normalize R appropriately. It is 
enslaved to take ri while it must play the role of an im- 
portant building block in the BLM pattern at the same 
time. Generating all b coefficients through Adler func- 
tion except bi and 62 would be a much more flexible way 
which reduces the infinite set of 02 observables to a set 
of few a2S. 



Every 02 observable connects to R through a specific 
channel. Temporarily we put the BLM pattern constraint 
aside. If we consider 02 to be the Adler function corre- 
sponding to i?, a very simple dependence of the 4th order 
coefficient of R to both 64 and 64 is obtained. So either 
one can predict both 64 and 64 or just predict 64 and use 
Adler transformation for 64 at the same time. Perform- 
ing the first choice for R results in a pattern crossing the 
exact pattern around Uf = 5 and for the second choice a 
crossing at n/ = 2. 

A more reasonable step would be to start with a spe- 
cific type of relation between two observables and then 
equip the situation with a BLM pattern. Of course it 
would be a waste of time if this ends to direct relations 
where we encounter observables whose predictions are 
equivalent. What we mean by a direct relation is a kind 
of linear relation between two variables. The situation is 
clarified later. For now we should start with a Crewther 
like relation between two observables R and 02: 

00 00 
[ro -I- riQs + ^ rifl/] [rg - riUs + ^ biUs'] 

i=2 i=2 

^C+^. (.) 

a 

where for the first observable we have R = riag + 
Si^2 ^^'^ the second one we have 02 = — t'iAs + 

Y^iL2^i^s^- Pa is the famous /3-function — X^i^o Z^**^'^^' 
and K ~ Y^^i^i^-^ ■ {^o---''i} and {Ki...Ki-i} deter- 
mines any coefficient bi in a special manner, considering 
all known coefficients prior to and Ki^i, 

1 1 C A 

bij = -rij - ——K^,_i)j + C^J i <i-l 
Tf 

bi[i-i)^ -''i(i-i) - -^K(i_i)(i_2)- (6) 

Despite the fact that R does not contain tq, it plays a 
crucial role in normalizing coefficients in Eq.®. To 
involve rp in the problem, one should perform prediction 
one order backward and normalize the result with the 
appropriate tq. Finally this will make ei coefficient be 
determined in terms of {ri...ri_i}. Taking K coefficients 
as the free variables, Eq.® reduces into the very simple 
direct relation bi = — -I- fiijif). The result comes out 
of such a relation suffers the same problem as a normal 
prediction for since bi and can exchange roles. The 
same is true about predicting ri through 6^. Although 
we are able to determine 6^ prior to using Eq.(l5]), at 
the same time 6^ has a direct relation with ri i.e., 6^ = 
CiTi -|- hiijif). So if anything is to give us real different 
results, it must have an indirect relation with r^: 

j-i 

bij— ''^[kj'Tij + mjjif(i_i)(j_i)] + dj, 
j=o 

i-l 

6i = J2^ii^^3 +hi{nf). 
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We have done a kind of converting such direct relations 
into indirect ones in the past when we tried to interpret 
DPA as changing weights of r^j coefficients in the normal 
BLM pattern. Obviously Eq.® is also a kind of indi- 
rect relation between and Ki^i- Following this line, it 



is clear that involving ith order K coefficients in Eq.® 
instead of {i — l)th order ones could be obtained by con- 
sidering in Eq.dSJ: 



[ro -I- ria + (r2o + r2i/3o)a^ + 



^ria''][ro-ria + {b2o + b2il3o)a^ + y^ha''] = C 



i=3 



i=3 



(7) 



where we've written all the coefficients in the /3- 

representation so n = r(i_i)(i_2) A-2 + J2]=o '^ijPo- The 
intrinsic property of this constraint would be to eliminate 
{Kq, Ki, K21, K32, K34} from the pattern. It is possible 
to use Eq.© alongside Eq.®. When get involved 
in Eq.©, an oversimple indirect relation appears. The 
indirectness comes from the absence of Ki parts in bio : 

bio = ~rio + Dij 

bij = -nj — - +Eij 0<j<i-l 

bi{i^i)= . (»j 

As in Eq.® we combined all lower order coefficients' 
dependences in constant parts Dij and Eij. Again pre- 
diction will depend on vq which should be resolved like 
before. 

One could rewrite Eq.® in terms of the coefficients 
of the /3-independent polynomials Vn in Eq.Q using the 
relations in [2ll| . The importance of doing this is that 
the idea in [21| strongly reveals existence of flavor- inde- 
pendent polynomials Vn- Conformal symmetry breaking 
totally happens through /3-function in this way. So it 
is more appropriate to rewrite Eq.([6]) in terms of these 
polynomial coefficients, even though this does not affect 
the directness of Eq. © . 



III. CONCLUSIONS 

Flexibility of any optimization prescription in pertur- 
bative analysis is the key to check the possibility of im- 
proving the prediction by taking deviation of lower or- 
der predictions into account i.e., finding patterns for the 
perturbative series. Taking DPA as finding parts of the 
observable that generates the appropriate BLM pattern 
for each order might help us in categorizing observables 
and their responses to the prediction respecting their de- 
viation patterns. But it could not be fully utilized to 
improve the prediction since for any order n we will suf- 
fer lack of information about {n — 1) parameters. In fact 
the effect of the symmetry group parameters on the BLM 
pattern is much more noticeable where the overall nor- 
malization plays a key role. A reasonable way to improve 
predicted patterns is to perform predictions based on in- 
direct relations. Observables are connected in a direct 
channel when they are expressed in terms of each other. 
In this case even specific insightful constraints such as 
the Crewther relation does not guarantee a simple indi- 
rect connection between observables. As a note for the 
reader, clearly we are not referring to the physical content 
of Crewther relation but just the mathematical form of it 
as a connection channel. So finding a mother constraint 
that produces elegant indirect relations is the challenging 
problem in the context of finding patterns. 
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